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We investigate pion superfluidity in the frame of two flavor Nambu-Jona-Lasinio model beyond 
mean field approximation. We calculate the thermodynamics to the next to leading order in an 
expansion in the inverse number of colors, including both quark and meson contributions at finite 
temperature and baryon and isospin density. Due to the meson fluctuations, the Sarma phase which 
exists at mean fleld level is washed away, and the Bose-Einstein condensation region at low isospin 
density is highly suppressed. 
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^ ■ I. INTRODUCTION 

' The study on Quantum Chromodynamics (QCD) phase structure is recently extended to finite isospin density 
. The physical motivation to study QCD at finite isospin density and the corresponding pion superfluidity is related 
' to the investigation of compact stars, isospin asymmetric nuclear matter and heavy ion collisions at intermediate 
OO , energies. 

■ While the perturbation theory of QCD can well describe the properties of new QCD phases at extremely high 

(-H ' temperature and density, the study on the phase structure at moderate temperature and density depends on lattice 
O QCD calculation and effective models with QCD symmetries. The lattice simulation at finite isospin chemical po- 
' I tential shows that there is a phase transition from normal phase to pion superfluidity phase at a critical isospin 
chemical potential which is about the pion mass in the vacuum. The QCD phase structure at finite isospin den- 
(— I sity is also investigated in low energy effective models, such as the Nambu-Jona-Lasinio (NJL) model [3] applied to 

I f'. quarks 0,i5|, 0,0, i] which is simple but enables us to see directly how the dynamic mechanism of isospin symmetry 

breaking operates. Near the phase transition point, the chiral and pion condensates calculated in this model are in 
[ good agreement with the lattice simulation 0. 

■ In a pion superfluid at zero baryon chemical potential, the quark and antiquark of a condensed pair have the same 
\^ ' isospin chemical potential and in turn the same Fermi surface. When a nonzero baryon chemical potential is turned on, 
it can be regarded as a Fermi surface mismatch between the quark and antiquark. The pion superfluidity in baryonic 
matter is recently discussed at mean field level in the NJL model in chiral limit [9*1 and in real case with finite current 
, quark mass [lO]. The pion superfluid can exist when the baryon density is not very high, otherwise the system will 
' be in normal phase without pion condensation because of the too strong mismatch. Inside the pion superfluid, the 
condensed state is separated into two phases. At small isospin chemical potential /xj, the homogeneous and isotropic 
Sarma phase (llj is free from the Sarma instability [ll[ and magnetic instability [12] due to the strong coupling and 
^ ■ large enough effective quark mass, it is therefore the stable ground state. At large /i/, while the Sarma instability 
can be cured via fixing baryon density ub to be nonzero, its mag netic instability implies that the inhomogeneous and 
anisotropic Larkin-Ovchinnikov-Fudde-Ferrell (LOFF) phase [ij] is favored than the Sarma phase. In the intermediate 
fj,i region, the stable ground state is the Sarma phase at higher ub and LOFF phase at lower ub- 

The Bose-Einstein condensation - Bardeen-Cooper-Shriffer (BEC-BCS) crossover at finite baryon and isospin chem- 
ical potentials is investigated in the NJL model [iJI ■ The pion condensation undergoes a BEC-BCS crossover when the 
isospin chemical potential increases. The point here is that the crossover is not triggered by increasing the strength of 
attractive interaction among quarks but driven by changing the isospin density. It is found that the chiral symmetry 
restoration at finite temperature and density plays an important role in the BEC-BCS crossover. 

Most of the work in the NJL model is mainly based on the mean field approximation to the quark mass and on 
the random phase approximation (RPA) for the Bethe-Salpeter equation for the meson masses hj. If one examines 
the thermodynamic potential in the mean field approximation, one sees immediately the deficit of this approach, 
viz, that only the quarks contribute to the thermodynamic potential with mesons playing no role whatsoever. This 
is clearly inadequate and unphysical, since one expects at least that the pionic degrees of freedom should dominate 
the system at low temperature, while the quark degrees of freedom should be relevant only in the chiral symmetry 
restoration phase. As such, this indicates that calculations in the NJL model must be performed beyond the mean 
field approximation. In Ref. the thermodynamics of a quark- meson plasma is calculated to order 1/Nc in an 
expansion in the inverse number of colors, and pions as Goldstone particles corresponding to spontaneous chiral 
symmetry restoration do control the thermodynamic functions at low temperature and density. 

A characteristic feature of the Sarma phase is the intermediate temperature superfiuidity (T7| : the superffuidity 
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happens at finite temperature but disappears at zero temperature. Since the mean field treatment is a good ap- 
proximation only at zero temperature [l8|, a careful study on the Sarma phase needs to go beyond the mean field. 
As for the BEC-BCS crossover induced by the change in density, the description on the BEC phase at low density 
should be closely related to whether the meson fluctuations are included or not. In this paper, we investigate the pion 
superfluidity in the frame of the NJL model beyond mean fleld approximation. We will focus on the effect of meson 
fluctuations on the Sarma phase and the BEC-BCS crossover at flnite temperature and baryon and isospin density. 

The paper is organized as follows. In Section |TT] we present the thermodynamics of the pion superfluidity and the 
gap equations for the chiral and pion condensates in the NJL model in and beyond mean fleld approximation. In 
Section Hm we calculate the phase diagram and see the meson effect on the Sarma phase and BEC-BCS crossover. We 
summarize and conclude in Section HVl 



II. THERMODYNAMICS OF THE PION SUPERFLUIDITY 



The two flavor SU{2) NJL Lagrangian density is defined as 

C = i^{ij''d^-mQ + n-fo)tp + G (VS^TstV')^ (1) 

with scalar and pseudoscalar interactions corresponding to a and tt excitations, where ip is the quark field, mp the 
current quark mass, G the coupling constant with dimension (GeV)^^, and ^ the quark chemical potential matrix in 
flavor space fi = diag{fiu, ^id) = diag{pB /S + fii /2, fiB /S — fxi /2) with fiB and /i/ being baryon and isospin chemical 
potential. The Lagrangian density has the symmetry [/^(l) S'C//(2) 5t/yi(2) corresponding to baryon number 
symmetry, isospin symmetry and chiral symmetry, respectively. However, at nonzero isospin chemical potential, 
the isospin symmetry SUi(2) breaks down to [//(I) global symmetry with the generator I3 which is related to the 
condensation of charged pions. At zero baryon chemical potential, the Fermi surfaces of u{d) and anti-c?(M) quarks 
coincide and hence the condensate of u and anti-d quarks is favored at sufficiently high /i/ > and the condensate of 
d and anti-u quarks is favored at sufficiently high /i/ < 0. We introduce the chiral condensate, 

(2) 

and the pion condensate, 

TT = V^i-tpij^T+ip) = V^irpi-f^T^ip) (3) 

with T± — (ri ± iT2) A nonzero condensate a means spontaneous chiral symmetry breaking, and a nonzero 

condensate tt means spontaneous isospin symmetry breaking. 

In mean field approximation the thermodynamic potential includes the condensation part and the quark part, 

n.mf = Gia^ + TT^) + flq, (4) 

and the quark part can be evaluated as a summation of four quasiparticle contributions [lo| . 
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where uji are the dispersions of the quasiparticles, 



UJ3 = E++ flB/i, UJi = E+ - flB/i (6) 



with the definitions 



E± = J{Ep ± /i//2)2 + 4G27r2, Ep = vV+^, (7) 



and the function g{x) is defined as g{x) — x/2 + rin(l -I- e~^/"^). The effective quark mass m is controlled by the 
chiral condensate, m — toq — 2Ga. The gap equations to determine the condensates a (or quark mass m) and tt can 
be obtained by the minimum of the thermodynamic potential f7m^(r, /i^, ///, m, tt), 

9^raf „ dn,nf d'^fljnf d'^ilmf , , 

om 07T om^ on'' 
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From the first order derivatives, we have 

1 dilq \ mo 



4G dm^ J AG 
.(G+ff)^0. (9) 

Considering the relations between /i^ and /ij, /i„ = /is/S + /i//2 and /i^^ = — /is/S + /i//2, the baryon and 
isospin density = —dilmj /dfiB and n/ = —dflmf /dfii can be expressed in terms of the u and d quark density 
n„ = -drimf/dfXu and = -dVlmf / d^i^, 

= ■^(»^n-«d)- (10) 

It is easy to see that ub plays the role of density asymmetry for pion condensation. For isospin symmetric matter 
with ub — Q the only possible homogeneous and isotropic pion condensed state is the BCS state. The Sarma state 
appears only in isospin asymmetric matter with ub 7^ 0. The two gap equations ^ and two number equations (jlOp 
determine self-consistently m, tt, /i/ and ^b as functions of T, n/ and nB at mean field level. 

We now consider the meson contribution to the thermodynamics of the system. The meson modes are re gard ed as 
quantum fluctuations above the mean field in the N JL model and can be calculated in the frame of RPA [15'| . For 
the mean field quark propagator with off-diagonal elements in flavor space, 

we must consider all possible channels in the bubble summation in RPA. In the pion superfluidity region, cr and 
charged pions are coupled to each other and the uncharged pion is decoupled from them. Using matrix notation for 
the meson polarization function 1 — 2Gn(g) 7], 
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with the quark bubbles 

^jk{q)^i J j^;^Tr[r*S{p + q)TkS{p)] , j,k^a, 7r+, 7r_, ttq (13) 
where the trace Tr = TrcTrpTrjj is taken in color, flavor and Dirac spaces and the meson vertexes are defined as 
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the meson masses are determined by 

det [1 - 2Gn(go + = Mj,q = 0)] = (15) 

with meson chemical potentials fi^ = 0, — M/i M;r_ = —fJ-i, fJ-Tia = 0- 

When the contribution from the meson fluctuations is taken into account, the total thermodynamic potential to 
order 1 /Nc in an expansion in the inverse number of colors becomes 

n^Qmf+^fl, (16) 

where the mean field part flmfjT, hb, fJ-ii tt) is shown in ([4]) and the meson part Qfi{T, fiB, fJ'ijmjTr) is expressed 
in terms of the polarization function TG], 



^fi = / tAi In det [1 - 2Gn(g)] . (17) 
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As is expected physically, the mesonic or collective degrees of freedom play a dominant role at low temperature, while 
the quark degrees of freedom are most relevant at high temperature [lal ■ 

In the chiral symmetry restoration phase at high temperature and/or high density, mesons are not stable bound 
states, but rather resonant states. They will decay into their quark-antiquark pairs. As a consequence, the determinant 
in the logarithm of (|17[) is a complex function in the meson energy plane and the imaginary part can be expressed 
as a scattering phase shift associated with quark-antiquark scattering. From the calculation in the NJL model with 
only chiral dynamics fl^, the meson width is small around the critical temperature but becomes remarkable when 
the meson mass is much larger than two times the quark mass, and correspondingly, the contribution from the phase 
shift to the thermodynamics is negligible at low temperature but significant when the temperature is high enough. 
For our calculation in the pion superfluidity phase, it can be estimated that the phase shift will be important in the 
BCS state at high density but its contribution is weakened in the BEC state at low density. Since we focus in this 
paper the Sarma phase and the BEC state which exist at low isospin density, we take pole approximation and neglect 
the scattering phase shift to simplify the numerical calculations. In pole approximation, the meson contribution can 
be greatly simplified as a summation of four quasiparticles, 

%i = ^r!„ (18) 

j 



Si, = 



with meson energies Ej — ^ Mj + . 

While mesons do not change the baryon density of the system, the charged pions modify the isospin density when 
the meson contribution to the thermodynamics is included, 

71/ = i (n„ + nj) + (n^^ - , (19) 

where Utt^ = —dflrr^/dfiTr^ and nT^_ — —dVl.„_/d^Tr_ are the 7r+ and 7r_ density. 

Up to this point, the order parameters m for chiral phase transition and tt for pion superfluidity have been regarded 
as the values minimizing fim/, and $7/; has been evaluated at these mean field values, m = rrimf and tt = nmf- While 
this is a correct perturbative expansion above the mean field, we may ask the questions: What is the feedback from 
the mesonic degrees of freedom to the order parameters and whether we could improve on these mean field values by 
regarding m and tt as variational parameters of the total thermodynamic potential Q = Qmf + f^/z? We now perform 
this procedure and see what the difference between the new and mean field condensates is. 

Taking the first order derivatives of the total thermodynamic potential with respect to the unknown quark mass m 
and pion condensate tt, we obtain the following modified gap equations, 

1 dflq dflfi\ Too 



AG dm^ dm? J AG 

In comparison with the mean field gap equations the fluctuation part in the thermodynamic potential leads 
to a new minimum at rumf+fi and nmf+fi that now differs from the mean field one at rrimf and nmf- It is easy to see 
that the structure of the new gap equations guarantees the two phase transitions. From the second gap equation for 
pion superfluidity, the trivial solution n — corresponds to normal quark matter, while the nonzero solution from the 
zero of the bracket corresponds to the energetically favoured pion condensed state. In the chiral limit, there are also 
two solutions of the first gap equation corresponding respectively to the chiral symmetry breaking and restoration 
phase. 

We now expand the fiuctuation part of the thermodynamic potential around the mean field minimum, 

nfl{T,tlB,^^I,m^,^T^)^ ^ ^ g(^2)» g(^2)j- ^/'(^' t^B, A^J, TO^ tt^) ^ ^ - ml^fY - tt^,/)" , (21) 



mf 



and, to further simplify the calculation, we consider the expansion only to the first order derivatives. Inserting the 
expansion into the new gap equations yields the following gap equations 

1 dflq \ Too 
icT ^d^J " 4G' 
dng 



G. + ^UO (22) 
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with two effective coupling constants Ga and G-^ defined by 

lid 
AG^^AG^ ^Q:^^f^'^^^^^B,^iI.ml,f,^Tl,f), 

^ mf 

d 

= G + ^-^^Jl{T,^MB,^J■I,■ml^f,'^l,f)■ (23) 

In comparing this group of coupled gap equations with the mean field one one observes that, in the chiral limit, the 
two groups take the same form, differing only in the effective coupling constants. The coupling constants in the scalar 
and pseudoscalar channels are the same at mean field level, but they become different and depend on temperature 
and charge densities when one goes beyond the mean field. If we take G„ = G,r = G, we recover the mean field case. 
That is, in this approach, the contribution from meson fluctuations is fully included in Ga and G^r. 

The above approach describes the thermodynamics of a quark-meson plasma with both chiral phase transition 
and pion superfluidity phase transition beyond the mean field at finite temperature and baryon and isospin density. 
The two new gap equations (j22p determine simultaneously the order parameters a (or to) and tt of the two phase 
transitions. In the chiral limit, the two phase transitions are fully separated from each other [3]: the chiral symmetry 
is automatically restored in the pion superfluidity phase. That is, the two order parameters do not coexist in the 
system. In the real word, chiral symmetry is not fully restored at any isospin chemical potential. However, a is much 
smaller than tt in the pion superfluidity region Q • Since we focus in this paper on the fluctuation effect on the pion 
superfluidity, we will, for the purpose of simplification in numerical calculations, neglect the a fluctuations and keep 
only the tt fluctuations in the gap equations. Namely, we take Ga = G in the following. 



III. PHASE DIAGRAMS IN AND BEYOND MEAN FIELD 



Since the NJL model is non-renormalizable, we should employ a regularization scheme to avoid the divergence in 
the gap equations. The simplest and normally used way is to introduce a hard three momentum cutoff |p| < A. In 
the following numerical calculations, we take the current quark mass toq = 5 MeV, the coupling constant G — 4.93 
GeV~^ and the cutoff A = 653 MeV This group of parameters ensures the pion mass TOjt = 138 MeV and the 
pion decay constant /^r = 93 MeV in the vacuum. 

In the treatment above, we considered the meson fluctuations as a perturbation around the mean field and took only 
the first order derivatives in the effective coupling constants ((23|) . If this treatment is good, the difference between the 
two pion condensates calculated in and beyond the mean field approximation should be small. To check the validity 
region of this method, we show in Fig[T]the two condensates as a function of isospin density at fixed temperature and 
baryon density. At low isospin density which corresponds to the BEG region, the difference between the two is really 
small, but it grows with increasing density and becomes large in the BCS region. Therefore, the approximation with 
only first order derivatives is good for the study of BEG, but the contribution from the higher order derivatives may 
be important for the BGS state. 
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FIG. 1: The pion condensates in mean field approximation (thin line) and including meson fluctuations (thick line) as functions 
of isospin density at fixed temperature and baryon density, go is the chiral condensate in vacuum and nj is the critical isospin 
density of pion superfluidity. 

The phase diagrams of pion superfiuidity in T — nj plane at fixed baryon density and in Ub — nj plane at fixed 
temperature are shown in Fig[21 The thin and thick solid lines are respectively phase transition lines in and beyond 
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mean field approximation which separate the normal quark matter at high temperature or high baryon density from 
the pion superfluidity matter at high isospin density. For pion superfluidity, the averaged Fermi surface of the paired 
quarks is controlled by isospin chemical potential and the mismatch is served by baryon chemical potential. Therefore, 
the Sarma phase which is induced by the Fermi surface mismatch may enter the pion superfluidity at nonzero baryon 
density. In mean field approximation, by analyzing the four quasiparticle dispersions w,;, the possible types of Sarma 
state and their thermodynamic and dynamic instabilities are discussed in detail in Ref . [7] . It is found that the Sarma 
phase is the ground state of the pion superfluidity at low isospin chemical potential. Very different from the BCS 
phase structure where the temperature of the pairing state is always lower than the temperature of the normal state, 
the Sarma phase appears in an intermediate temperature region and the normal state exists in lower and higher 
temperature regions, see the mean field phase transition line in the T — nj plane at low isospin density in Fig|21 
However, the phase structure in the T — nj plane is significantly modified when the meson fluctuations are included. 
From Fig[21 the meson effect reduces greatly the pion superfluidity region, and the critical temperature is suppressed 
from about 150 MeV in mean field treatment to about 80 MeV in the case beyond the mean field. A qualitative 
change resulted from the meson fluctuations is that the intermediate temperature superfluidity or the Sarma state in 
mean field calculations is totally washed away, and the normal quark matter is always above the BCS pairing state. 
In the ns — ni plane, the phase diagram in mean field approximation is similar to the one in ub — fJ-i plane obtained 
in Ref. pLOj . and again the meson effect reduces remarkably the pion superfluidity region. 




rii/no rii/no 

FIG. 2: Tlie pliase diagrams of pion superfluidity in T — nj plane at fixed baryon density ns/no = 0.5(left panel) and in 
riB — ni plane at fixed temperature T — 50 MeV(right panel) in mean field approximation (thin lines) and including meson 
fluctuations (thick lines), no = 0.17//m^ is the normal nuclear density. The solid lines are the phase transition lines and the 
dashed lines are the BEC-BCS crossover lines. 

It has been argued both in effective theory and lattice simulation that at finite but not very large isospin density 
and zero baryon density, the QCD matter is a pure meson matter, i.e., a Bose-Einstein condensate of charged pions. 
At ultrahigh isospin density, the matter turns to be a Fermi liquid with quark-antiquark cooper pairing[l|. Therefore, 
there should be a BEC to BCS crossover when the isospin chemical potential increases. There are some equivalent 
quantities to describe the BEC-BCS crossover induced by changing charge density pJl]. Among them are the root- 
mean-square radius of the Cooper pair which is small in BEC and large in BCS, the s-wave scattering length which is 
positive in BEC and negative in BCS, the condensate scaled by the Fermi energy which is large in BEC and small in 
BCS, and the fermion chemical potential which is negative in BEC and positive in BCS. In the following we take the 
chemical potential to characterize the BEC-BCS crossover. For relativistic pion superfluidity, the chemical potential 
which controls the BEC-BCS crossover is ^//2 — m depending on temperature and baryon density through the 
effective quark mass m, and fii can be viewed as the binding energy of the bound state of quark and antiquark in the 
BEC limit. In Fig|3]we show /i//2 — m as a function of nj at flxed temperature and baryon density in and beyond 
mean field approximation. In both cases, the effective chemical potential goes up from negative to positive values 
with increasing isospin density. The zero point, namely the BEC-BCS crossover point, is located at nj/no = 1.68 
in mean field treatment and nj/uQ — 1.29 in the case with meson fluctuations. The crossover lines determined by 
fij /2 — m = in T — nj and ub — nj planes are shown in Figl^l When the mesonic fluctuations are included, not 
only the pion superfluidity region is greatly reduced, but also the BEC region is strongly shrunk. 

In the BCS limit of the pion superfluidity, the isospin density is high and the paired quark and antiquark is 
weakly coupled. At the critical temperature, the condensate disappears and the weakly coupled fermions are excited 
separately, and the system is a Fermi liquid. In the BEC limit, however, the isospin density is low and the paired 
quark and antiquark is tightly coupled. In this case, above the critical temperature, the system becomes a Bose 
liquid of tightly bound pions, and the quarks should be too heavy to be excited. This means that, at the critical 
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FIG. 3: The effective chemical potential /ii/2 — m as a function of isospin density at fixed temperature T = 50 MeV and 
baryon density nslna — 0.5 in mean field approximation (thin line) and including mesonic fiuctuations (thick line), no is the 
normal nuclear density. 



temperature mesons are lighter than quarks in the BEC limit and quarks are lighter than mesons in the BCS hmit. 
To confirm the BEC-BCS crossover picture obtained above by calculating the effective chemical potential inside the 
pion superfluidity, we show in FigH] the meson mass and quark mass m as functions of isospin density at the 

critical temperature and fixed baryon density. With increasing isospin density, the quark mass drops down but the 
meson mass goes up monotonously. The two lines cross at about n//no = 1.4 which is qualitatively in agreement 
with the BEC-BCS crossover value determined by /z//2 — m = 0. 
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FIG. 4: The meson mass Af^^ (solid line) and quark mass m (dashed line) as functions of isospin density at the critical 
temperature Tc and fixed baryon density ns/rzo = 0.5. no is the normal nuclear density. 



IV. SUMMARY 



We have investigated the thermodynamics of a pion superfluid at finite isospin density in the frame of two flavor 
NJL model beyond the mean field approximation. Considering the fact that mesons, in particular pions because of 
their low mass, dominate the thermodynamics of a quark-hadron system at low temperature, the mesonic fluctuations 
should be significant for the phase structure of pion superfluidity. By recalculating the minimum of the thermodynamic 
potential including meson contribution, we derived a new gap equation for the pion condensate which is similar to 
the mean field form but with a medium dependent coupling constant. From our numerical calculations, the main 
effects of the meson fluctuations on the phase structure are: 1) the critical temperature of pion superfluidity is highly 
suppressed and the Sarma phase which exists at low isospin chemical potential in mean field approximation is fully 
washed away, and 2) the BEC region at low isospin density is significantly shrunk. 
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